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Spectral and temporal coincidence and single-particle photon wave packets are described and their 
widths and durations are found. The degree of entanglement is characterized by the experimentally 
measurable parameter R defined as the ratio of the coincidence and single-particle spectral widths. 
In the frequency representation, this parameter is found as a faction of the duration r of the pump 
pulses. The function R(t) is shown to have a minimum and even in the minimum, at rather natural 
conditions, the parameter R is found to be very high J? m i n ~ 73 3> 1. The Schmidt number K 
is found for both short and long pump pulses and interpolated for arbitrary pulse durations. All 
functional dependences of R and K are found to be identical and numerical difference between 
them is shown to be not exceeding 20%. Two-time temporal wave function of a biphoton state is 
investigated in details, and a rather significant difference between the cases of short and long pump 
pulses is found to occur. In the case of long pulses, the temporal parameter R t is defined as the 
ratio of durations of the single-particle and coincidence signals, and Rt is shown to be very close to 
the Schmidt number K. 
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I. INTRODUCTION 



Entanglement is a fundamental feature of multipartite 
systems characterizing the degree of correlations between 
particles. Entangled states find their wide application 
in quantum information, quantum cryptography and re- 
lated topics. Spontaneous Parametric Down Conversion 
(SPDC) is one of the most often used processes for pro- 
duction of entangled biphoton states. The most popu- 
lar and widely discussed type of entanglement of SPDC 
states is the entanglement with respect to polarization 
variables. On the other hand, in addition to polariza- 
tion, SPDC photons are characterized also by other de- 
grees of freedom, such as angles determining directions of 
propagation and lengths of wave vectors (or frequencies) 
of photons. In contrast to polarization, these degrees 
of freedom belong to the class of continuous variables. 
Multi- and bipartite states with continuous variables are 
quite attractive because they correspond to an infinitely 
high dimensionality of the Hilbert space and, potentially, 
their degree of entanglement can be very high. But in 
practice, of course, there are limitations. Nevertheless, 
there are conditions when the degree of entanglement 
with respect to continuous variables can be really high, 
much higher than the maximal achievable degree of en- 
tanglement with respect to discrete variables. For ex- 
ample, in Refs. [l[ SPDC biphoton states were shown 
to be highly entangled with respect to angular variables. 
In this work we investigate entanglement of SPDC bipho- 
tons with respect to their frequency variables, or spectral 
entanglement. 
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We characterize the degree of entanglement by the re- 
cently introduced Q parameter R defined as the ratio of 
widths of single-particle and coincidence wave packets. 
We find R analytically in two asymptotic cases of short 
and long pump pulses (with the parameter separating 
these two regions explicitly defined) and we show that 
in these two cases R is extremely high (~ 300 or more). 
By using both interpolation formulas and exact numer- 
ical calculations we show that in the case of the type-I 
phase matching we consider, R is very high (~ 70) even 
in the most unfavorable case of intermediately long pump 
pulses (~ 1 ps) . 

In contrast to all other entanglement quantifiers, the 
parameter R is rather easy experimentally measurable. 
For double-Gaussian bipartite states the parameter R is 
known Q to coincide with such a well known entangle- 
ment quantifier as the Schmidt number K . However, the 
wave functions characterizing the biphoton frequency dis- 
tributions are rather far from the double-Gaussian ones. 
Nevertheless, as we show in this paper, even for such 
functions the parameters R and K are very close to each 
other. By using specific features of the biphoton wave 
functions we calculate in Section V the Schmidt number 
explicitly in both cases of short and long pump pulses and 
we find that in both cases all functional dependences of R 
and K are identical, and there is only a small difference 
between them in numerical coefficients. This means that 
the parameter R is, indeed, an appropriate quantifier of 
entanglement which can be both found theoretically and 
measured experimentally. 

In addition to spectral analysis we describe also a tem- 
poral structure of the SPDC biphoton wave packet deter- 
mined by its two-time wave function. This structure has 
many very interesting features. We derive a series of very 
simple formulas, which are used to explain and character- 
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ized analytically parameters of the calculated temporal 
envelopes of single-particle and coincidence signals. We 
find that features of these curves are qualitatively differ- 
ent in the cases of short and long pump pulses. In the 
case of short pump pulses duration of coincidence tempo- 
ral signals is shown to be varying in very large limits with 
varying time of observation, and this makes the temporal 
picture with short pump pulses inappropriate for defining 
the temporal parameter R, R t . In contrast, in the case of 
long pump pulses, duration of coincidence temporal sig- 
nals is shown to be independent of the observation time, 
the parameter R t is well defined and is shown to be very 
close to the Schmidt number K . 



II. GENERAL FORMULAE 

SPDC is the process in which a classical pump wave 
propagates trough a birefringent crystal, where some 
photons of the pump are absorbed and give rise to the 
birth of two photons of smaller frequencies, signal and 
idler. If ki and k2 are wave vectors of emitted photons, 
the QED state arising in the SPDC process is given by 
the sum of the vacuum state and superposition of two- 
photon states 
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where |ki, k 2 ) = a k a k | vacuum), a k is the operator of 
creation of a photon in a state with the momentum (wave 
vector) k, and ^(ki, k 2 ) is the biphoton wave function in 
the momentum representation. General expressions for 
this wave function are well known [1, Q , and for pump 
pulses of a finite duration the biphoton wave function can 
be written in the form 

*(kij_, k 2 ±; ui, lo 2 ) =X J dr J dt J dk p ± J duj p 

xE p (kp±, top) exp[i(k p _L - ku_ - k 2 _i_) • rj_ + iA ■ z] 
x exp[i(wi + uj 2 - w p )t], (2) 

where \ is the crystal susceptibility, indices p and _L in- 
dicate, correspondingly, the pump and the plane perpen- 
dicular to the laser axis (z-axis), E p is the Fourier trans- 
form of the electric field amplitude of the pump, and A 
is the longitudinal mismatch 



2 2 



k 2 - 



1^1 _j.2 _ 

A 1-L 



1^2 i,2 

„2 K 2-L' 

(3) 

n p , m, and n 2 are the refractive indices of the pump and 
emitted photons. 

Let us make now a series of simplifying assumptions 
and approximations. Let us consider the type-I SPDC 
decay process e — > o+o, in which both of emitted photons 
belong to the ordinary wave and, hence, n\ = n 2 = n Q . 
Integration over time t in Eq. @ from — oo to oo gives 



2n5(ujp — u>i — u> 2 ), with which the integral over oj p can 
be easily taken. Similarly, in the wide-crystal approxi- 
mation the limits of integration over can be extended 
to — oo and oo in both transverse directions which gives 
(27r) 2 (S(k p j^ ki^ k2j_), and this can be used to take 
the integral over k p j_. Let us restrict our consideration 
in this work only by the case of a purely collincar prop- 
agation of emitted photons, in which ki j_ = h 2 j_ = 
and, hence, k p ^ = 0. Besides, let us consider here only 
the case of degenerate central frequencies cj^ and w 2 °\ 
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0J2 = where uiq is the central frequency 

of the pump. We will assume also that deviations from 
these central frequencies, = uj\ y2 
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are small: 



|^i,2| ^ wo- Finally, we will take the spectral amplitude 
of the pump in a Gaussian form 
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where r is the duration of pump pulses and lu p = u>i + 
u 2 — (Jo + v\ + v 2 . 

Under the formulated assumptions and approxima- 
tions Eq. ^ for the biphoton wave function is reduced 
to a much simpler form 

*(wi,wa)«£i,(a;i+wa) [ dze iA ^' u *> (5) 
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or, after integration over z, 

^(wi, lu 2 ) oc E p (uji + lo 2 ) sine 
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where sinc(u) = sinu/it, L is the length of the crystal 
(along the laser axis), and the phase mismatch A(u>i, u 2 ) 
is determined by the equation following from Eq. ([3]) 



A = k p — k\ — k 2 
{u)\ + Lu 2 )n p (uJi + lu 2 ) — u)\rio{u)\) — u)2no(u) 2 ) 
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Note, that as we will be interested mainly in shapes and 
widths of photon spectral distributions, in Eqs. ([5]) and 
([6]) and henceforth we drop all the constant coefficients 
on the right-hand sides of such equations by using the 
proportionality symbol. If not specifies differently, all 
curves describing spectral or temporal distributions will 
be normalized by the condition of being equal to unity 
at maxima. 

Owing to the assumption about small deviations of 
frequencies u>i, 2 from lu^\, the expression on the right- 
hand side of Eq. ([7]) can be expanded in powers of Pi and 
v 2 with only two lowest orders to be taken into account 
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where A is the dimcnsionlcss constant characterizing the 
temporal walk-off 
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are the group velocities of the pump and or- 



dinary waves (v g °^ > vf ) ; B in Eq. ||SJ) is the dispersion 
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For LHO3 crystal and Ao = 400 nm numerical values of 
the constants A and B are: A = 0.17 and B = 0.069. 
Though small, nevertheless, the quadratic term in Eq. 
([5]) is crucially important for determining finite width 
of the photon single-particle distribution (see the next 
section, part B). Note also that in Eq. ([5]) we have 
dropped the dispersion arising from the wave vector of 
the pump and proportional to k'{(y\ + v 2 ) 2 . This term 
is dropped because it depends on the same combina- 
tion of frequencies as the linear term v\ + v 2 and de- 
termines only a small correction to the linear term not 
changing qualitatively its functional dependence on v\ 
and v 2 . In the same approximation we can further sim- 
plify the dispersion term taken into account in Eq. |8|): 
v\ + 2/f = + V 2 f\ + (l/! - v 2 f\ -> \{v x - v 2 f. 

With this substitution the final expression for the bipho- 
ton wave function takes the form 
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III. COINCIDENCE AND SINGLE-PARTICLE 
DISTRIBUTIONS (SHORT PULSES) 

A. Coincidence spectrum 

The coincidence spectrum is determined by the 
squared wave function v 2 )\ 2 of Eq. (TiT]) at a given 

value of v 2 . For practical purposes, it is more convenient 
to consider the coincidence and single-particle probabil- 
ity densities as functions of wavelengths Ai,2 rather than 
frequencies v\ t % = 2nc/Xi t2 — ttc/Xq. For the crystal 
length and pump-pulse duration equal to L = 0.5 cm and 
t = 50 fs, the coincidence spectral distribution dw^°' /dXi 
is characterized by narrow curves of Fig. [T] The wide 
curves in the pictures of Fig. [T] describe the pump spec- 
trum. The curves of the pictures (a) and (b) are calcu- 
lated at the second-photon wavelength A2 taken equal to 
800 nm and 870 nm. 

Fig. [1] clearly shows that (i) the coincidence spec- 
trum is much narrower than that of the pump and 



FIG. 1: Coincidence and pump spectra (narrow and wide curves) 
as functions of Ai at A2 = (a) 800 nm and (b) 870 nm. 



{ii) the coincidence spectral width AA^ does not de- 
pend on A2. If we take the wavelength of the second 
photon exactly equal to the double pump wavelength, 
A2 = 2Ao = 800 nm, the peaks of the pump and coinci- 
dence spectra coincide and are located at Ai = 800 nm 
(Fig. [T]Ja))). If we shift A 2 to the right, both pump 
and coincidence spectra move to the left from the point 
A2 = 800 nm. However, for the coincidence spectral dis- 
tribution this shift is slightly larger than for the pump. 
For this reason the amplitude of the coincidence signal 
appears to be somewhat smaller than in the case of lo- 
cation at 800 nm. However, the width of the coincidence 
spectral distribution remains practically the same as in 
the case A2 = 800 nm. Hence, the coincidence spectral 
width AX\ is a constant parameter, independent of tun- 
ing of the second-photon detector. 

Analytically the coincidence spectral width can be eas- 
ily found directly from Eq. (fTl"]) . In the case of short 
pulses the Gaussian function in this equation is much 
narrower than the sinc-function. For bipartite wave func- 
tions having the form of a product of two functions, one 
of which is much narrower than another, the coincidence 
distribution is determined by the narrower factor, i.e. in 
our case by the sinc-function. Evidently, it's width with 
respect to varying v\ at a given value of v 2 is given by 
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where the factor 2.78 is the full width at half-maximum 
(FWHM) of the function sinc 2 (u). Numerically, at L = 
0.5 cm, in terms of wavelengths, both analytical expres- 
sion and calculations shown in Fig. [1] correspond to 

(c) 

AX\> = 0.658 nm. For comparison, FWHM of the pump 
spectral function cxp[— {vi + z/ 2 ) 2 T 2 /41n2] at a given v 2 
equals to 
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which gives (in terms of wavelengths, at r = 50 fs and 

A = 400 nm) AA^ p) = 2A ! Ai/ lpunip /7rc = 18.8 nm. In 
other words, the pump spectrum considered as a function 
of Ai, 
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is 28.57 times wider than the coincidence spectrum. 

As shown below, the ratio of the "sine" and "pump" 
widths [([T2"]l and fT5]) ] determines the parameter separat- 
ing the regions of short and long pulses 
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The last expression on the right-hand side of this equa- 
tion shows that 77 equals the ratio of the double pump- 
pulse duration to the difference of times during which a 
crystal is crossed by the pump and idler/signal photons. 
Also r] can be considered as the dimcnsionless pump-pulse 
duration. In SPDC the pump pulses are short if 77 -C 1 
and long if t] ^> 1. The boundary pulse duration oc- 
curs at r\ ~ 1. At L = 0.5 cm the equality r\ ~ 1 yields 
t = 1.43 ps. As shown below, r\ is also the key control 
parameter, determining a value of the entanglement pa- 
rameter R. 

Non-coincidence of the peaks of the pump and coinci- 
dence spectra seen in Fig. [T](6) deserves a special com- 
ment. For clarification of its origin, it's worth to return 
to the frequency pictures. In Fig. [5] the coincidence and 
pump spectra are shown as functions of i^i under the 
conditions when the detector registering the second pho- 
ton of the pair ( "2" ) is tuned to frequencies shifted as 
compared to wq (or v 2 shifted from zero) (a) to the blue 
and (b) to the red sides. Both the pump and coincidence 
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FIG. 2: Coincidence and pump spectra dw^ (u\) / dv\ at (a) V2 = 2 
and (b) U2 = —2, v\ and vi are in units of 10 fs. 

spectra also appear to be shifted from luq (v\ = 0), but 
in the directions opposite to the shift of lo 2 . As is seen, 
the peaks of the pump and coincidence spectra do not 
coincide at sufficiently large values of these shifts. Inter- 
esting enough, in both cases the coincidence peak appears 
to be located at the right wing of the pump, which corre- 
sponds to the left wing in the wave-length picture of Fig. 
[TJ This lack of symmetry between the pictures (a) and 
(6) of Fig. [2] is related to the dispersion and, as we will 
see in the following section, it is related to asymmetry of 
the single-particle spectrum. Mathematically the results 
discussed here follow directly from the definition of the 
phase-mismatch A, e.g., taken in the form ([5]). Location 
of the coincidence spectral peaks is determined by the 



condition A = 0. This gives a quadratic equation of 
solution of which is given by 
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The first and second terms on the right-hand side of this 
equation have different parity. They are, correspond- 
ingly, odd and even functions of u 2 , and this explains 
the lack symmetry between the pictures (a) and (b) of 

Fig.rj 

The approximation used in the expansion of the square 
root in Eq. (|15p consists in the assumption 1 1 <C 
Ajo/8-B, which is satisfied for parameters we use in this 
paper. Actually, if here the approximate expression for 
^1(^2) was needed only for getting a simple analytical 
formula explaining the exact calculations of the spectra 
in Figs. [T]and[2l later this will be shown to be crucially 
important, in particular, for calculations of the Schmidt 
number (Section V). 



B. Single- particle spectrum and the entanglement 
parameter for short pump pulses 

The single-particle photon spectrum is determined by 
the squared wave function of Eq. (jlip integrated over v 2 



dw( s > 
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As the exact analytical calculation of the integral over 
v 2 is hardly possible, let us perform this integration ap- 
proximately with the narrow sinc 2 -function substituted 
by the 5-function &\v — i>(i>i)], where v = v\ + v 2 and 
v{v\) is given by equation similar to Eq. (|15p 
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By using now v as the integration variable in Eq. (|16| 
(instead of v 2 ), we find immediately the following expres- 
sion for the single-particle spectrum of emitted photons 
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In Fig. [3] the spectral distribution determined by Eq. 
(|18[) is shown together with the pump spectrum (corre- 
spondingly, the wide and narrow curves). The FWHM 
of the wide curve in Fig. [3] equals to 
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FIG. 3: Single-particle (wide) and pump (narrow) spectra (densi- 
ties of probability of the corresponding photon distributions). 



Numerically, at r = 50 fs, Eq. JTSJ) yields AA^ } = 
195 nm, which is 10 time larger than the pump spectral 
width. 

The parameter R defined as the ratio of the single- 
particle to coincidence spectral widths is easily found 
from Eqs. (|12|) and (TT9|) to be given by 
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where the subscript "short" emphasizes that this result is 
valid only for short pump pulses, 77 <§C 1 , and the control 
parameter rj is defined in Eq. (|14| . The numerical coef- 
ficient 55.06 in the last expression of Eq. (|2"D"|) was found 
at the same values of all parameters which were used 
earlier, except the pump-pulse duration r which was not 
fixed yet. At r = 50 fs (?] = 0.0348) Eq. flU} yields 
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This degree of entanglement is really high. As it's seen 
from the given derivation, such an extremely high de- 
gree of entanglement is related to two accompanying ef- 
fects: narrowing of the coincidence and broadening of 
the single-particle spectra of SPRC photons (compared 
to the pump spectrum). The second of these two effects 
is a specific feature of the type-I SPDC process. Indeed, 
in the case of the type-II pase matching (o — > e + o), 
in the linear approximation, the expansion of the phase 
mismatch in powers of v±, takes the form 
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As in the case of the type-II phase matching the emit- 
ted photons "1" and "2" are different (o- and e- waves), 
derivatives of their wave vectors over frequencies are dif- 
ferent too, k[ k 2 . For this reason, the phase mismatch 
of Eq. (|22p already in the linear approximation appears 



to be depending on both the sum and difference of fre- 
quencies, vi + v 2 and v^ — v 2 . In this case the linear ap- 
proximation in u\, v 2 is sufficient for determining finite 
widths of both coincidence and single-particle distribu- 
tions. If k[ — k 2 is of the same order as k[ + k' 2 , the 
difference between the single-particle widths and that of 
the pump spectrum is not too large. I.e., the effect of 
broadening of the single-particle distribution is missing 
or is only weakly pronounced. In contrast to this, in the 
case of the type-I phase matching both emitted photons 
belong to the same ordinary type of a wave and, hence, 
fcj = k' 2 . In the linear approximation the dependence of 
the phase mismatch A (|22|) on the difference of frequen- 
cies v\ — v 2 disappears, and the linear approximation ap- 
pears to be insufficient for determining the finite width 
of the single-particle distribution. To define this width 
one has to take into account the second order (disper- 
sion) term in the expansion of A in powers of V\, v 2 . As 
this term is rather small, the width of the single-particle 
distribution appears to be rather large. I.e., equality of 
k[ and k' 2 gives rise to the above-described broadening 
of the single-particle spectrum compared to that of the 
pump. This effect provides at least the 10-fold increase 
of the parameter R compared to what can be achieved in 
the type-II SPDC process. 

Note that the above-described difference between the 
type-I and type-II SPDC processes was understood 
rather long ago Q). New elements of our considera- 
tion consist in demonstration of importance of this ef- 
fect for producing byphoton states with extremely high 
degree of entanglement and in explicit evaluation of the 
achievable degree of entanglement. As far as we know, 
there was only one experiment where both coincidence 
and single-particle SPDC spectra were measured simul- 
taneously [2j- However, the difference between the coin- 
cidence and single-particle distributions observed in this 
work was not too high and corresponded to the ratio of 
widths on the order of 2. We assume that one of the rea- 
sons why this ratio was so small is in the type of phase- 
matching (type-II). 

In Ref . H the degree of spectral entanglement for the 
type-II SPDC process was evaluated with the help of de- 
composition of the wave function in a series of products 
of Schmidt modes, which is different from our approach. 

In some works (see, e.g., Ref. @) the 2D photon distri- 
bution was measured and plotted in the (wi, 0J2) plane. 
However, as far as we know, such results were never used 
to evaluate quantitatively the degree of entanglement. 

Finally, in this work we discuss only possibilities of di- 
rect measurements of the coincidence and single-particle 
spectral and temporal distributions of photons. We do 
not discuss here connections between results obtained in 
such a way with interference measurements, e.g., in Rcfs. 
0, HO, El El El- We hope to return to this problem 
elsewhere. 
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IV. LONG PUMP PULSES AND 
INTERPOLATION FOR ARBITRARY PULSE 
DURATIONS 

Until now we have assumed that the sine-function in 
Eq. (jlip is much narrower than the Gaussian function. 
This is absolutely true, e.g., for the pulse duration r = 
50 fs, which was chosen for estimates and illustrations. 
However, if the pulse duration is significantly longer, the 
relation between the two functions on the right-hand side 
of Eq. may be reversed: the Gaussian function can 
become much narrower than the sine-function. In this 
case the coincidence distribution is determined by the 
pump spectrum. E.g., for i/ 2 = it has the form 
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In the integral (|16|) determining the single-particle proba- 
bility density the narrow pump spectrum can be approx- 
imated by 8(i>2 + v-i) to give 
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The coincidence and single-particle width found from 
Eqs. (H3J and (J2J]) are equal to 
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and the corresponding entanglement parameter R is 
given by 
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The entanglement parameter i?i on g (|26| is a linearly 
growing function of the pulse duration r. E.g., at r = 7ps 
Eq. ([H yields r) w 4.87 and i? long w 269 > 1. 

For evaluating the entanglement parameter in the 
whole range of variation of the pump-pulse duration r 
(from short- to long-pulse regions) we use the simplest 
interpolation formula 



R(r) 



(27) 



Note that there is an amusing, though approximate, co- 
incidence relation between numerical coefficients in Eqs. 
([20)1 and J2BJ: V2tt In 2/2.78 w 1.004 V2/78^/(2 5 / 2 In 2). 
Hence, with a very good accuracy we can take coefficients 
in Eqs. ([2H| and ([2"6")1 coinciding with each other to get 
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FIG. 4: The parameter R and K vs. the control parameter r) (1141 1 . 



The function i?(?7) is plotted in Fig. [4) The curve has 
a minimum at 770 = 2 -1 / 3 = 0.79, which corresponds to 
t = 1.14 ps, and in the minimum 
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So, the parameter i?(r) is very large and the degree 
of entanglement is very high in both earlier considered 
cases of short and long pump pulses, r <C tq and r ^> 
to. But R(t) is large and entanglement is high also in 
the most unfavorable case of intermediate pump-pulse 
durations around To. Actually, the value of R(t) around 
its minimum is determined by the well defined parameter 
equal to the square root of the ratio " length of the crystal 
divided by the pump wavelength" . This result shows 
that in the type-I degenerate collinear SPDC process, 
practically, spectral entanglement is never small. 



V. THE SCHMIDT NUMBER 

As mentioned above (in Introduction), in the case of 
double-Gaussian wave functions, the parameter R coin- 
cides identically with the Schmidt number K Q, which 
is widely recognized to be a good quantifier of the degree 
of entanglement in pure bipartite states. However, the 
wave function of the form ([IT)) is not double-Gaussian. 
Actually, if the pump spectral amplitude can be taken 
Gaussian, the sine-function with the argument containing 
both linear and quadratic terms in variables 1^2 cannot 
be even modeled by any Gaussian form. For this rea- 
son, the question that often arises is whether the param- 
eter R can be used as the entanglement quantifier for 
states characterized by such non-double-Gaussian wave 
functions, while relation between R and K remains un- 
known. The question is reasonable and, to be answered, 
it requires the Schmidt number to be explicitly calcu- 
lated. Solution of this problem is given below. 

Let us use the well known general integral definition 
of the Schmidt number for bipartite systems with con- 
tinuous variables [Eq. (10) of Rcf. c 3|]- In the case of 
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spectral entanglement this definition takes the form 
K = N 2 J dvxdv<xdv' x dv'<£Sl{yx, v-z)^l{y x , v' 2 ) 

where N is the norm of the wave function 

2 



(30) 



N = dv x dv 2 ll'Oi, v 2 )X 



(31) 



and ^(^i, v 2 ) is given by Eq. (fTTjl . 

In a general case, the integrals in Eqs. ([30]) . (f3Tj) can- 
not be calculated analytically, and even their numerical 
calculation is a rather difficult problem. To perform in- 
tegrations we will use approximations following from the 
main features of the wave function pip described above. 
Let us consider separately the cases of short and long 
pulses, ?)Cl and 77 S> 1. 

1) Short pump pulses. 

The mismatch A(^i, v 2 ) (|8|) entering the argument 
of the sinc-function in the definition of v 2 ) is 

quadratic in v\. Hence, the equation A = has two 
solutions, Vi (v 2 ) and v^{y 2 ). Only one of these solu- 
tions is compatible with the condition \v\. 2 \ <C and 
this solution is given by equation similar to Eq. (JT5J) 



v x {v 2 ) = f{v 2 ) = -v 2 + 



AB 
Auq 



In the vicinity of v x the mismatch A(^i, v 2 ) can be ap- 
proximated by the linear function of v\ 

A = A(i>i, v 2 ) w — 1 — /(1/2)] ■ (32) 
c 

The same approximation can be applied to all other func- 
tions ^> in the definitions of the Schmidt number (|30|) and 
norm of the wave function pip , which take then the form 



x exp 

x exp 

x exp 

x exp 
and 



K~ 


^N 2 = 


{vi 


+ V2?T 2 




31n2 


(^1 


+ v' 2 ) 2 t 2 




3 In 2 


K 


+ vifr 2 



sine 



81n2 



;ln2 



smc 



dv\dv 2 dv' x dv' 2 

vi - Mi) 

v'x - M)) 

v'x-M<x) 



N 



dv\dv 2 exp 



LA 
~2c~ 
LA 
~2c~ 
LA 
~2c~ 
LA 
~2c 



(1/1 + v 2 ) 2 t 2 



(33) 



41n2 



x < smc 



LA 
~2c 



V\ - /(^2)J 



(34) 



Still, the integrals ([33)1 and (|34[) are too complicated for 
analytical calculation. But we can remember now that 
in the case of short pulses (cr <C AL/2) the exponential 
functions in all integrals are much wider than the sinc- 
functions. Hence, with a good accuracy, we can take off 
the exponential functions from the integrals at values of 
their arguments determined by the condition that the 
corresponding phase mismatch turns zero. With such a 
procedure applied to Eq. ([34]) we get 

45W 



27TC . 

Is = -j—^ I dv 2 exp 



A 2 uj 2 In 2 



2 3/: 



— (ln2)V4 r 



^0 
ABt' 



(35) 



where we took into account that sinc 2 (x)dx = ir and 

JZo e ~ X dx = 2r ( 5 / 4 ), r denotes the gamma- function, 
and r(5/4) = 0.9064. 

For calculation of integrals in Eq. (|33|) we use the 
following exact equation for the integrated product of 
two sinc-functions 



dxsmc(x) sinc(x + y) = 7rsinc(y). (36) 



Known or not, this equality can be easily proved, e.g., by 
means of analytical continuation to the complex plane x 
and integration by the residue method. With the help of 
the rule (|36p . wc perform integration in Eq. (|33[) sepa- 
rately in the first and second pairs of lines under the sym- 
bol of integrals (correspondingly, over vi and with 
the exponential functions taken out of the integrals, as 
explained above. As a result we get 



K^N 2 



f2ircY 

[la) 



dv 2 dv' 2 exp 



A 2 uj 2 In 2 



smc 



LA 

^ V2 - V2 > 



(37) 



Note that in the last expression the difference vi(v 2 ) — 
v\{v 2 ) in the argument of the sinc-function is approxi- 
mated by v' 2 — v 2l i-c, here the quadratic term of the ex- 
pression in Eq, (fT5j) is dropped both in f(v 2 ) and f{v 2 ). 
This is possible because in Eq. (|37|) the argument of 
the exponential function differs from that of the sinc- 
function, and hence, small quadratic terms in the sinc- 
function argument are not needed for making integrals 
converging, as this was, e.g., in the case of integrations 
over v\ and v>[. The integrals in Eq. (|37p are easily 
calculated with the help of the same procedure as used 
above, with the slow exponential function taken off from 
the integral over v' 2 at v' 2 ~v 2 . The results are given by 



K~ l N 2 



_ (2nc 
= \LA 



(2-kc 
\LA 

3 



dv 2 exp 



(2 In 2 



A 2 lu 2 In 2 



1/4 



Au 
Bt 



(38) 
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and 



K. 



N 2 



short 



(21n2) 1 /4 /5\ A 3 / 2 L 



K^ 1 N 2 
= 0.785- 



.4 



47 v / crAo 
57.5 



(39) 



Comparison with Eq. (|2H)) shows that in the case of 
short pump pulses all functional dependences of i? s hort 
and A'short are identical and, numerically, the difference 
between i? s hort and -Kshort is very small, ~ 4.5%. 
2) Long pulses. 

In the case of long pump pulses (77 3> 1) calculation of 
the Schmidt number is similar to that described above for 
short pulses, though the roles of the exponential (pump) 
and sine- functions are reversed: the pump spectral func- 
tion is narrow and the sinc-function is wide. Owing to 
this, the sinc-function in Eq. (jllj) can be simplified in a 
way different from that used in Eqs. (|32|). (|33|) and (|34|) . 
In the case of long pulses the sum V\ + v 2 is very small 
and the linear term in the argument of the sinc-function 
can be dropped. As the result we can write down Eqs. 
(|30]) and (J3T]) in the form 



x exp 
x exp 
x exp 
x exp 



K- l N 2 = J dvidv 2 dv' 1 di/ 2 

'LB {v x -v 2 f 



iyx - 


f V2?T 2 
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31n2 


{v\ - 




\ 


3 In 2 


K- 


f V 2 ) 2 T 2 



smc 



smc 



81n2 
8 In 2 



smc 



2c 




LB { Vl 




2c 


Uq 


LB (v[ 


~V2? 



smc 



2c oj 
LB K 



2c 



and 



N 



dv\dv 2 exp 



wo 



+ V 2 ) 2 T 2 



(40) 



4 In 2 



x < smc 



LB (i>i - v 2 ) 2 



2c 



(41) 



The norm N is calculated easily with v 2 = —V\ substi- 
tuted into the argument of the wide sinc 2 -function, after 
which the integral over v 2 in (|4ip takes the Gaussian form 
to give 



2\/^Iri2 . 
N = — / dv x 



(42) 



(ILBvl 
c 

V CUJQ 
2\J TT In 2 / CLOq 4y / 7T 

2LB~3~' 



where the last factor on the right-hand side of Eq. (|42|) 
arises from integration over the dimensionless variable 
x = u i \/2LB/c u>o : dx sinc 2 (x 2 ) = 4y / 7r/3. 



The product of two exponential functions in the first 
and second lines under the symbol of integrals in Eq. ([40|) 
equals to 



exp 



!ln2 



[{vi + v 2 f + {vi + v 2 ) 2 



exp 



4 In 2 



V 2 + Vn 



(y2-v 2 ) 2 



(43) 



In dependence on v\ this narrow function has a maximum 
at v\ — —(y 2 + i/' 2 )/2, and this is the value that has to be 
substituted instead of V\ into the wide sinc-functions in 
the first and second lines under the symbol of integrals in 
Eq. (|40p . This substitution makes the sinc-function inde- 
pendent of v\ , and integration of the exponential function 
(|43p over v\ is easily performed. Integration over i/[ in 
the last two lines of Eq. (j40|) is performed exactly in the 
same way. Altogether, these two integrations (over v\ 
and v[) reduce Eq. ([4^]) to a simpler form 



r ^ 1 »r9 47rln2 , 
K~ 1 N 2 = — / dv 2 dv' 2 exp 



Sln2 



{v2-v' 2 y 



smc 



LB 

8c wo 



(3^2 + vtf 



smc 



LB 

8c wo 



{v2+?>v' 2 y 



(44) 



where the first and second sinc 2 -functions arise, corre- 
spondingly, from integration of the first and second pairs 
of lines in Eq. (|4"U1) . In Eq. (144j) . again, the exponential 
function is narrow and sinc-functions are wide in their 
dependence on, e.g., v' 2 . Owing to this, we can substi- 
tute v' 2 by v 2 in the arguments of the sinc-functions and 
take easily the arising Gaussian integral over v' 2 to get 



r _ 1n2 = 8V2(n In 2) 3 / 2 



dv 2 sine 



A f2LBi% 
\ cw 



8V2(7rln2) 3 / 2 [^Tq 



2LB 



^i (2 3/ 2 _ j 
105 v 



(45) 



where the last factor in square brackets is the exact ex- 
pression for the integral (ixsinc (x 2 ). 

Eqs. (|42|) and (|45|) determine the Schmidt number in 
the long-pulse regime 



105^ 



72 V21n2 (2 3 / 2 - 1) y/BLXo 
Arj FT 
^B\Xo 



0.6 



44?y. 



(46) 



Comparison of this result with that of Eq. (|26|) shows 
that in the case of long pulses, as well as in the case 
of short pulses [Eqs. ([20]) and ([39]) ]. the Schmidt num- 
ber and the parameter R are rather close to each other. 
Again, all functional dependences of Eq. ([46]) and (|26|) 
are identical. The parameters Ai ong and i?i on g differ only 
by numerical coefficients. In the case of long pulses this 
difference is ~ 20%, i.e., somewhat larger than in the 
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case of short pulses but, still, small enough to consider 
R as a good entanglement quantifier. 

The Schmidt number for arbitrary pulse durations can 
be evaluated with the help of the same square-root inter- 
polation as in the case of R [Eq. ([28)) ] : 




(47) 



The function K(rj) is plotted in Fig. [4] (the dashed curve) 
together with R(ij). The difference between these two 
curves is more pronounced in the case of long pump 
pulses and is almost negligible in the case of short pulses. 
But even in the case of long pulses the parameters R and 
K are close enough to use the experimentally measurable 
parameter R as the entanglement quantifier. 

The ratio of R and K parameters (|2"5)l and (|4T|) equals 

to 



Kjrj) 



1.04 




(48) 



Numerical coefficients on the right-hand side of this equa- 
tion are not related to any parameters of the medium or 
pump and arise form In 2, ^pK and similar factors. For 
this reason the ratio K/R can be considered as a univer- 
sal function depending only on the control parameter rj 
and valid for any crystals (with the type-I phase match- 
ing) and any pump wave lengths. Moreover, Eq. (|48[1 and 
all other results derived above [except the very final nu- 
merical estimates in Eqs. (gOl), (H3) , (ED- ©, 63), and 
(|47p ] can be applied to any processes generating bipartite 
states of the form (fTT|) with arbitrary pairs of variables 
substituting v\ and v 2 and with any values and physical 
origin of the constants A and B. In particular, the results 
described here can be applied directly to the earlier con- 
sidered angular entanglement of SPDC biphoton states 
in a continuous pump laser p| . Specifically, for establish- 
ing one-to-one correspondence between the spectral and 
angular entanglement, we have to substitute frequencies 
^1.2 by angles 9\^ 2 determining directions of the photon 
wave vectors kj. 2- I n accordance with the approach of 
Refs. we can consider the angular divergence of the 
pump ao as a varying parameter substituting in the case 
of a monochromatic pump the pump-pulse duration r. 
Then, to reduce the angular wave function in the "par- 
allel geometry" of Refs. pj to the form of Eq. we 
have to substitute the constants A ((9|) and B (fT0|) and 
the control parameter 77 (|14|) by 



A= — ^ 
Ao 



A 



41n2A 



2n p Ao ' 



V = 



n p 



TiaoL n' 



(49) 



where now n' p is the angular derivative of the pump re- 
fractive index n p . With these substitutions done, we can 



use all the results derived above for the Schmidt number 
K and the parameter R. In particular, with with the 
help of Eq. (f2"5)l we can find the minimal value of the 
parameter R characterizing angular entanglement in the 
parallel geometry at varying angular divergence of the 
pump ao: 



R 



angu 
min 



lar 



2L 
j p Aq 



(50) 



This result was not obtained earlier in refs p|, and it 
shows that the degree of angular entanglement is high 
even in the most unfavorable circumstances owing to the 
same large factor as in the case of spectral entanglement, 

1. 

Finally, Eq. (|48|) and the picture of Fig. [4] can be used 
for direct experimental determining of the Schmidt num- 
ber (if the 20%-accuracy evaluation of the degree of en- 
tanglement provided by R is insufficient). To find K one 
has to measure R, to calculate the value of the control 
parameter rj (|14|) corresponding to experimental condi- 
tions, to find its location at the plot of Fig. U] and to 
find K . Alternatively, with known R and rj the Schmidt 
number is easily calculated with the help of Eq. 



VI. TEMPORAL STRUCTURE OF BIPHOTON 
WAVE PACKETS 

A. General formulae 

Temporal and coordinate features of biphoton states 
are characterized by the function ^(t\, t 2 ] Zi, z 2 ) de- 
pending on two times t\ and t 2 and two coordinates, z\ 
and z 2 , 



t 2 ; zi, Z2) = J dw\dw 2 L02) 
x cxp j— % u)\t\ - k\{u)\)z\ + uj 2 t2 - k 2 {uj2)z2 I • (51) 

The wave function ^ determines the probability density 

t 2 ; zx, z 2 )\ 2 (52) 



dw 



dt\dz\dt 2 dz2 

of finding (registering) photons "1" and "2" in small 
vicinities of arbitrary points [z\, t±) and (z 2 , t 2 ) in the 
coordinate-time plane (z, t) under the condition that z\ 
and z 2 are within the crystal, i.e., < zi 2 < L. 

The exponential factors depending on (zi, t\) and 
(z 2 , t 2 ) on the right-hand side of Eq. (|5l"j) are the time- 
coordinate wave functions of one-photon states with mo- 
menta fci and k 2 (or frequencies u)\ and uj 2 ) 



i(a>2*2— £222) _ 



(z 2 , i 2 |fci) = (z 2 , t 2 \a\ \vac) 



(53) 
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The two-time and two-coordinate wave function of Eq. 
(foTj) is the projection of the two-photon state vector given 
by the second term on the right-hand side of Eq. |T]) upon 
the product of two coordinate states 

t 2 ; 21, z 2 ) = 
(zi, h\(z 2 , t 2 \ *(ki, k a )|ki, k 2 ) (54) 

ki,k 2 

for the case ki )2 _l = 0. 

Note that owing to the presence of coordinate parts 
in these single-photon wave functions, transition from 
spectral to temporal picture determined by Eq. (|5ip 
is not equivalent to simple double Fourier transform of 
^(lui, uj 2 ) used in Refs. @, [HI]. Though some results of 
these two types of calculations can coincide (see Eq. 
and explanations below), in a general case transition to 
the temporal picture is inseparable from transition to the 
coordinate representation, and this is the transformation 
determined by Eq. (|51[) that gives a rigorously defined 
two-time and two-coordinate biphoton wave function. 

As it was done previously, the wave vectors ki^i^i, 2) 
in Eq. (f5"Tj) can be expanded in powers of small differ- 
ences v\ t 2 — — 2~' JJ o^ with only the first and second 
orders to be retained 



fci, 2 (wi 



^0 
2c 



Vl,2 
J°) 



2B 

CLU 



(55) 



By using this expansion, the definition of Eq. (|5ip . and 
the expression of Eq. ([5]) for the frequency-dependent 
wave function ^>(u>i, uj 2 ), we get the following integral 
representation for 

f ' ^ ) 2 ^~ 2 

$(ti, t 2 \ 21, z 2 ) oc dz d^d^cxp 1 - ^ 



x exp 



x exp 



.Az , , 
+ 1 — [yi + v-2 ) 
c 



. B / zi +z 2 
1 2 — [y x - v 2 ) 

CUlQ 



(56) 



where, as previously, the squared frequencies cf 2 = 
-Ayi — v 2 ± [y\ + v 2 )\ 2 in the expansion (|55[) are approxi- 
mated by \{y\ — v 2 ) 2 with the terms ±i(i/i — j^H^i +^2) 
and jft'i + v% ) 2 dropped as giving only small corrections 
to other terms oc v\ + f 2 and oc (1/1 + v 2 ) 2 in the expo- 
nential functions of Eq. (f56|) . 

In variables u = z^i + v 2 and v = V\ — v 2 the double 
integral over v\ and i/ 2 in Eq. (|56[) splits for the product 
of two Gaussian integrals over u and v which are easily 



taken to give 



*(ti, i 2 ; 21, 22) oc 



min(zi,22) ^ z 



2 In 2 



x exp 



exp 



Az 



c(ti + t 2 ) c(zi + z 2 ) 



(°) 



4BA (2i + 2 2 - 22) 



ti-t 2 - 



2v 



z\ - 2 2 



■1; 



(o) 



(57) 



Here we restricted the region of integration over 2 by 
^(21 + 22) (instead of L) because we interpret 2 as the 
coordinate of birth of a single-idler pair, and Z\ and 22 as 
possible locations of the emitted photons. As these pho- 
tons propagate to the direction of growing 2i )2 , these co- 
ordinates cannot be smaller than the coordinate of their 
birth, i.e., 2 must be smaller than min{zi, 22}. 

At the exit from the crystal, where z\ = z 2 = L, Eq. 
(l57l) takes the form 



*(*!, i 2 ) = *(ti, L; t 2 , L) cc 



fir 



Vl^ 



2 In 2 



x exp 



2 (L-z) (ii+t 2 ) 



-1 2 



,(e) 



,(°) 



X exp < I 



-*a) s 



8BA (£-2) 
Finally, by introducing new time variables 



(58) 



*i,2=ti, 3 t+ = ^(Ti+t 2 ), t-=h-t 2 (59) 

we can slightly simplify the first exponential factor on 
the right-hand side of Eq. j[5"8"|) and get 



*2) oc 



exp 



21n2 



(Az - ct + ) 2 ^ exp 



7TC t 



2 ,2 



8BX (L - 2) 



(60) 



The time shift L/vg"^ in Eq. (|59p is the time during which 
the signal and idler photons moving exactly with their 
group velocity cross all the crystal from to L. Hence, 
if £1 = t% = is the time when the peak of the temporal 
envelope of the pump enters the crystal, the zero point 
for tx,2 is the time when the signal and idler photons, 

emitted at ti = i 2 = and moving with the speed v g ° , 
reach the end of the crystal. 

Note that substitution of the integration variable 2 = 
z' + L and some additional shift of the time variables t\ 
and t 2 reduce Eq. (f60|) to the same form as that of Eq. 
(22) in the paper by Keller and Rubin This equation 
of Ref. @ was obtained with the help of a simple dou- 
ble Fourier transformation rather than via the transition 
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from momentum to time-coordinate representation de- 
termined by Eq. (|51[) . Though seemingly surprising, this 
coincidence is quite understandable. Indeed, in terms of 



-L and 0, and at the exit 
0. At this point, the one- 



z', the crystal borders are 
from the crystal z[ = z' 2 = 
photon wave functions of Eq. (|53|) (with z substituted 
by z') turn into e - ^ 1 - 2 * 1 - 2 , and the transformation of 
Eq. (|51[) appears to be equivalent to the double Fourier 
transformation. However, this is a very exceptional case 
related to the assumed in Ref. Q location of the crystal 
at (— L,0). Our attempt [lfl to apply a simple dou- 
ble Fourier transformation to the case when the crystal 
location is taken as (0, L) gave the result qualitatively 
different from that of Eq. (|60[) and, evidently, wrong. 
Here in Eq. (|60[) we correct it with the help of deriva- 
tion based on the transformation from the frequency to 
time-coordinate representation of Eq. (|51|) . 

Actually, the results of calculations have to be invari- 
ant with respect to the choice of coordinates determining 
location of the crystal. We can easily find that the trans- 
formation determined by Eq. (|5ip provides such invari- 
ance whereas the simple double Fourier transformation 
does not. Indeed, the difference between Eq. ([50]) and 
Eq. (10) of Ref. [3 is in a different location of the 
singularity in the integral over z: at z = L (the exit 
edge of a crystal) in Eq. (16"01) and z = (the entrance 
edge) in Eq. (10) of Ref. [15| . Location of the singu- 
larity is determined by the factor in front of (vi — v 2 ) 2 
in the last line of Eq. (|56[) and in similar equations of 
Refs. @ and [l5[. Let in a general case the coordinates 
of the crystal edges be z = L\ and z = L 2 = L\ + L. 
We are interested in the biphoton wave function at the 
exit edge, i.e., z± — z 2 ~ L 2 . In this case, the coeffi- 
cient in front of iy\ — v 2 ) 2 in the last line of Eq. ([55]) 
is proportional to z — L 2 , and the singularity in the fi- 
nal integral over z [similar to that of Eq. (f60|) ] occurs 
at z = L 2l i.e., always at the exit edge of the crystal. 
In contrast to this, in the scheme of calculation with the 
double-Fourier transformation the coefficient in front of 
(v\ — v 2 ) 2 in the last line of Eq. ([56]) becomes propor- 
tional to z, i.e., the singularity in the final integral over 
z occurs at z = 0. This point can be arbitrarily located 
with respect to the crystal and, evidently, the results are 
not invariant with respect to the choice of the crystal co- 
ordinates. The only case when the results of calculations 
in these two approaches coincide is the case L\ = —L 
and L 2 = 0, and this is just the case of Ref. [f|. 



B. Short pump pulses 

The probability density of registering photons at 
the exit from the crystal at times t\ and ti is given 
by the squared absolute value of the wave function 
^f(t\, t 2 ) (|60|) . determining a surface in the 3D space 
jti, t 2 , l^l 2 !- For short pump pulses r = 50fs and for 
the same values of all other parameters which were used 



above, the surface \P(ti, t 2 ) is shown in Fig. [5] 



|*£ti, ta) 
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FIG. 5: The function t 2 )\ 2 for r = 50fs, (a) 3D plot and 

(b) the top view. 

The surface |^(ii, t 2 )\ 2 is seen to have the form of a 
long plateau symmetric in t\, t 2 . The plateau is rather 
wide in the beginning (at small values £1,2) and is nar- 
rowing and heightening with a growing value of t\ + t 2 . 
Around the point t\ = t 2 = 285.25 the plateau turns 
into a rather high and sharp peak. The diagonal squared 
wave function |^(ti, ii)| 2 characterizes evolution of the 
plateau-peak height in its dependence on t\ + t 2 , and this 
function is shown in Fig. [6] In accordance with the 3D 
picture of Fig. E](a), the curve of Fig. [5] switches on 
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FIG. 6: The diagonal squared wave function ii)| 2 - 

rather quickly at t = (in a time on the order of the 
pulse duration of the pump r), and then it has a long 
and slowly growing plateau ending by a high and narrow 
peak at t = 2.8525 ps followed by a very quick turn-off. 

The temporal shape and parameters of coincidence 
and single-particle signals are defined in the same way 
as in the cases of frequency or angular distributions. 
The coincidence signal is determined by the probabil- 
ity density dw^ c \t\) / dt\ of registering at the exit sur- 
face of the crystal a photon "1" at a varying time t\ 
under the condition that the photon "2" of the same 
pair is registered at the same place at some given time 
t 2 , dw^it^/dh oc |$(ii, t 2 )\ 2 at t 2 = const. For 50-fs 
pulses of the pump, two examples of coincidence curves 
are shown in Fig. [7J These curves correspond to the 
very beginning of the plateau and to the peak of the sur- 
face t 2 ) | 2 in Fig. El The FWHM of the curves (a) 
and (b) in Fig. [7] are equal to At[ c) (t 2 = 0) = 486.3 fs 
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FIG. 7: Coincidence signals at (a) t 2 = and (b) t 2 = 2.8525, ps. 

and At^ (t 2 = 2.8525 ps) = 6fs. This huge difference in 
FWHM shows cleariy that in the temporal picture the 
coincidence widths At^ (or the duration of the coinci- 
dence signal) depends very strongly on the registration 
time t 2 of the second photon of a pair. This contrasts 
with the spectral picture where A^} was shown to be 
almost independent of v 2 . A strong dependence of At± 
on t 2 makes the temporal picture inappropriate for defin- 
ing the parameter R and for determining in such a way 
the degree of entanglement. Such a procedure can be 
used only in the spectral picture. Some analogies can be 
found in the coordinate and momentum pictures for mas- 
sive particles, where spreading of wave packets in the co- 
ordinate representation makes the latter inconvenient for 
evaluation of the degree of entanglement via the width- 
ratio parameter R. Whereas the momentum represen- 
tation is free from such a problem because momentum 
distributions do not spread 0- In the case of photons 
dispersion or diffraction divergency plays the role analo- 
gous to spreading. 

The single-particle temporal density of probability to 
register a photon by a single detector at a varying time t\ 
is determined by the integrated squared two-time wave 
function of Eq. 

^^oc J dt2 m u t 2 )\ 2 . (6i) 

For the pump-pulse duration as short as r = 50 fs the 
result of numerical integration is plotted in Fig. [5] Dura- 
tion of the single-particle signal, defined as FWHM of the 
curve in Fig. [5J equals to 2.837 ps, which is more than 50 
times longer than the pump pulse. Two other interesting 



dw^ /dt ± 




50 100 150 200 250 300 
ti [10 fs] 



FIG. 8: The normalized time-dependent single-particle probability 
density. 



features of the curve are its constant value 

at a very long interval of time and the anomalously long 
front wing (at small t\). Explanations are given below 
in the subsection D on the basis of simplified analytical 
formulas. 



C. Long pump pulses 

All the results shown in Figs. [5][8j correspond to a short 
duration of the pump pulses (i) C 1). In the case of long 
pulses (77 3> 1) the picture changes rather significantly. 
For t = 2 ps the 3D pictures of Fig [5] are substituted by 
those of Fig. [5J The surface \&(ti, t 2 ) is seen to have 




200 _ 400 C00 -200 200 _ 400 000 

t a U 



FIG. 9: The function |*(ti, t 2 )| 2 , 3D plot for r = 2ps. 

the form of a narrow long bump concentrated along the 
diagonal t\ = t 2 in the plane (t\, t 2 ). As seen clearly 
in Fig. [10](&), the width of the layer \&(ti, t 2 ) is almost 
constant in the middle part of the bump, and the layer 
narrows slightly at its ends. Quantitative characteristics 
of the temporal photon distribution will be derived and 
presented in the last part of the next Section. 



D. Analytical characterization 

a) Simplified formulae. If the pump pulses are short 
enough, cr <C AL, the Eq. the first exponential factor 
in the integrand of Eq. (|60[> ) is very narrow. Location 
of its peak determines the coordinate zq = ct+/A, a 
small vicinity of which gives the main contribution to 
the integral over z. Under these conditions we can sim- 
plify the integrand by substituting z = zq into the pre- 
exponential factor 1 / ^/L — z and by expanding l/(L — z) 
in the second exponential factor in powers of \z — zq\/zq: 
1/(L - z) w 1/(L - zq) + (z- z )/(L - z ) 2 . As a re- 
sult of such simplifications, the integral in Eq. (foT)]) takes 
the Gaussian form with finite limits of integration, which 
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gives rise to the sum of two error functions 



^/LA - ct + 



exp 



LA - ct A 



Erf 



-Erf 



/— — -AL - rf+ 
V21n2 - 



i a 



ct. 



ct- 



AL - ct + 

2 



where 



V21n2— +ia , 

r V LA - ct A 



ttctA 



16V21n2BA n 



(62) 



(63) 



The validity criterion of the approximation used for 
derivation of Eq. (f6"2"]) is \z— zq\ <C L—zq with zq = ct + /A 
and \z — zq\ ~ ct/A, which gives 



LA - ct. > cr. 



(64) 



As t + is associated with the location zq of the peak of the 
pump envelope in the crystal, the condition (|64[) means 
that the pump pulse is located far enough from the end 
of the crystal, at a distance exceeding significantly ct/A. 
Note, that the times t+ = and t+ = LA/c are asso- 
ciated with the pump pulse entry to and exit from the 
crystal. 

If, in addition to ([M]) . the time t+ or coordinate zq 
obey the conditions t + > t or zo > ct/A, the pump 
envelope is located far from both edges of the crystal. In 
this case the formula of Eq. (f6"5)) can be further signifi- 
cantly simplified. Indeed, if a very short pump pulse is 
located deep in the crystal, the limits of integration over 
z can be extended to — oo and +oo. This is equivalent to 
approximation of the first and second error functions in 
Eq. ([B"2"]) by +1 and -1, and reduces finally Eq. (|6"2"j) to 
the simplest form: 



1 



\JLA - ct+ 



exp 



ct. 



LA - ctj 



(65) 



Eqs. (f6"2"]> and (f6"S"|) are much more convenient than the ex- 
act formula (|60p for obtaining analytical expressions for 
parameters characterizing coincidence and single-particle 
signals. 

b) Localization of the temporal biphoton wave 
packet. The localization region of the wave function 
fy(tx, t 2 ) m the (t±, t 2 ) plane can be defined in the fol- 
lowing way: we plot graphs of |^| 2 as a function of t_ 
at all given values of f+ and find coordinates tt(t+) of 
these curves at half-maxima. Then the functions tt(t+) 
together with the lines t+ — and and t + ss 3ps deter- 
mine the borders of the |*| 2 localization region. With 
the functions ft (t + ) calculated numerically from the Eq. 
(|6T)|) the result is given by the diagram of Fig. [TQl in 




FIG. 10: Localization region of the function ^(ti, *2); 
1 1 , t%, 1 4. , and t — are in units of 10 — 14 s; (b) is the increased section 
III of (a). 



which t^t(f + ) and tZ(t+) determine the upper and lower 
borders of the shaded area. Defined in such a way, the 
diagram of Fig. agrees perfectly with the view-from-the- 
top of the surface t 2 ) in Fig. [5](6). In Fig. LTOlonc 

can distinguish clearly three regions of t+, I, II, and III, 
in which the dependence t_(i+) is significantly different. 
Below these regions are analyzed separately. 

c) Coincidence width in the region II. The most 
extended of three regions in Fig. [10] is the region of 
intermediate values of £ + , II, where the wave function 
^f(t+, t-) can be satisfactorily approximated by the sim- 
plest exponential formula of Eq. (|65[) . This equation 
gives 



±- 



4 /ln2ABA 



it Act 



(LA - ct+) . (66) 



By definition. Eq. (|66p determines the upper and lower 
boundaries of the localization region in Fig. [10] at a given 
value of t + . However, the same equation can be used to 
find the upper and lower boundaries of the localization 
region tf at a given t<i. To find ^(^2) one has to make a 
substitution in Eq. (|66|) t + = \(ti+t 2 ), to replace ft by 
t\—t%, and to solve the arising equation with respect to ti . 
In the first order in a small parameter 2ylri 2 BXq/ttct w 
0.1 C 1 (at t = 50 fs) the solution is given by 



tf(t 2 ) = h it 4 



ln2L?A 



it Act j c V ttct 



From this equation we obtain the following analytical 
expression for the width of the coincidence temporal dis- 
tribution of photons in the region II of Fig. [TU] 



Ati c) (t 2 ) = t+-«r = - 

c 



BX In 2 
ttAct 



(LA - ct 2 ) . (68) 



This expression shows that spread of photons coming in 
pairs to the exit from crystal over time of their arrival is 
determined by a joint influence of the temporal walk-off 
and dispersion. The spread is large at small values of 
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the observation time characterized by t 2 and falls with a 
growing t 2 . 

d) The region I. In the region I the simplest approx- 
imate formula of Eq. ([6"5]) is invalid. But the formula of 
Eq. ([62]) works perfectly and can be used to find ana- 
lytical expressions determining both the shape of bound- 
aries of localization region in Fig. 1 101 and the coincidence 
width of the temporal distribution of photons. As in the 
region I the time t+ is very small, the first error function 
in Eq. ([6"2"]) equals unity, and the term t + can be dropped 
everywhere except the first term in the argument of the 
second error function to give 



* cx e- y [1 + Erf(6 + zy 2 )], 



(69) 



where b = -\/21n2i+T, y = ct-y/a/LA, and a is defined 
in Eq. (|6"5|) . From Eq. (|6"9"|) we find easily (e.g., with 
the help of simple calculation in "Mathcmatica" ) that at 
b^= the HWHN (half-width at the half-maximum) of 
1^1 2 equals one, and at growing but small values of b it 
falls as 1 — 0.24 b. In terms of t+ and i_ this yields 



£(*+) = ± 



LA 



1 - 0.24 



V2 In 2 £, 



4£(21n2) 1 / 4 VABA^ / ^ y/lhrttj 



(70) 



The lines (t + ) determined by Eq. (JTUJ) are sharper than 
the lines, determined by Eq. I66[) . The crossing points of 



these lines occur at t+ ~ 43 fs and they determine the 
boundary between the zones I and II in Fig. [TO] As seen 
well in Fig. [TQ1 the width of the localization region in the 
£_-direction is maximal at t + = where, in accordance 
with Eq. (fTTTj) . it equals to 



At_ 



t+(t+=Q)-tZ(t + =Q) 



8L(21n2) 1 / 4 v / ABA^ 



1.46 ps. 



(71) 



jAi_ max , where At_ max is given 



e) Pulse shape of a single-particle signal. A large 
width of the wave-packet localization region in the t_- 
direction at t + = explains the origin and parameters 
of the anomalously long front wing of the single-particle 
signal shown in Fig. [S] Indeed, owing to the condi- 
tion t+ = 0, in the (t±, t 2 ) frame the points C a and Cp 
in Fig. [10] have coordinates ti a = — t 2 .u 
and hp = - i 2 , p : 

by Eq. ([7T]) . As the ti-dependent single-particle signal 
can be considered as the sum of all coincidence signals 
at all values of t 2 , it's clear that at t 2 < t 2a there are 
no signals at all, either in coincidence or single-particle 
measurements. Hence, t 2a is that tim since which the 
single-particle signal become detectable. Actually, t 2a 
is the beginning of the transient period during which the 
amplitude of the single-particle grows. This transient pe- 
riod ends at t 2 = t 2 p, and its duration equals to Ai_ max 
(I7ip . The same value characterizes the duration Tf w of 



the front wing of the single particle signal considered as 
a function of t\ : 



Tfw = tic — t\p = 



8L(21n2) 1 / 4 VABA^ 



(72) 



Another interesting feature of the single-particle signal 
is its constant value during all the its duration after end- 
ing of the initial transient period (Fig. [S]). This feature 
is easily explained by Eq. ([65]) and ([68]) . Indeed, at a 
given ti the single-particle yield ([6T]) can be estimated 
approximately as the product of the hight ti)| 2 
and width At 2 (ti) of \^ft±, t 2 )\ 2 considered as a func- 
tion t 2 . The first of these two factors is given by Eq. ([65]) 
with t- = and, evidently, h)\ 2 cx 1/{LA - cti). 

As for the width At^ c) (ti), it is given by Eq. ([68]) with 
t 2 on its right-hand side substituted by t\ , and this gives 
At 2 c \ti) oc LA — ct\. Hence, with a growing ti the peak 
hight of the function ^(ii, ^2 ) | 2 (in its dependence on 
t 2 ) grows, whereas its width falls, but their product re- 
mains constant (independent of t±), i.e., the dependence 
dw^ (ti) / dti has the form of a plateau. 

Finally, the total duration of the single-particle sig- 
nal is determined by the condition that the real part of 
the argument of the second error-function in Eq. 
changes its sign, i.e. 



Tsingle ~ n - (e) (n y (JO J 



L 



This result is absolutely understandable qualitatively: 
^single is the period of time between the arrival to exit 
from the crystal of almost first and almost last observ- 
able SPDC photons or the delay time between the arrival 
of the first ordinary-wave photons and of the pump (ex- 
traordinary wave) . 

f ) The region III in the diagram of Fig. 1101 This 
is the region behind the peak of the the photon distribu- 
tion in Fig. [5](a), where the amount of photons is vary 
small. Nevertheless, this region is interesting, because 
here formation of the localization region of the temporal 
wave packet is related to a mechanism qualitatively dif- 
ferent from those described above for the regions I and 
II. In the region III the condition of Eq. I|64p is not ful- 
filled, the fraction 1/(L — z) in the integrand of Eq. ([ST)]) 
cannot be expanded in powers of (z — Zq)/{L — Zo), and 
the simple formulas of Eqs. (|62[) and ([B3]) are invalid. 
But some estimates can be done directly on the ground 
of Eq. ([6"D]) . Indeed, from the second exponential factor 
in the integrand on the right-hand side of Eq. ([6T?]) wc 
can deduce the following estimate of the characteristic 
half-width of the wave-packet localization region 



At_± ~ ± 



^BX (L - z) 



(74) 



If the peak of the laser pulse is located exactly at the end 
of the crystal, ct + = LA, then the maximal deviation of 



15 



z from L, where the integrand of Eq. ([60]) is not small, 
is L — z ~ ct/A. Substitution of this expression into Eq. 
([74]) gives 



At: 



± 



VBXot 
VcA 



(75) 



This expression corresponds to the very beginning of the 
zone III in Fig. [TU1 

If ct + > L, the peak of the pump pulse is located 
behind the crystal and on a small part if its rear wing 
remains in the crystal. By expanding the Gaussian ex- 
ponent in powers of L — z, we find that now L — z ~ 
(?t 2 /\A(ct + — LA)]. By substituting this expression into 
Eq. (|74|) . we find the law by which the zone III in Fig. 
1101 slowly narrows 



At_ 



BXn 



A(ct+ - L) ' 



(76) 



g) Long pulses. The regime of long pulses of the 
pump arises if the pulse duration is so long that t 3> 
LA/c or 77 S> 1 where r\ is given by Eq. In this 

case the term proportional to z can be dropped in the 
argument of the first exponential function in Eq. (|60|) to 
give 



t 2 ) oc exp 



F(t. 



where 



F(t-) 



dz 



exp 



8BX (L - z) 



(77) 



(78) 



Eq. (|77[) shows that for long pulses the dependences of 
the wave function ^ on t + and t- arc factorizcd, and 
the dependence of \& on t + repeats the time-dependence 
of the pump envelope. As for the integral over z in Eq. 
([75)) determining the dependence of \& on it can be 
calculated to expressed F{tJ) in terms of the error func- 
tion. To show this, we have to substitute the variable of 



integration z = L 







1 - i 

s 







where 



To 



1 8BXpL 

C V 7T 



(79) 



and s is a new purely imaginary variable. Then the inte- 
gral (|78|) takes the form 



( 



F(t_) oc M 



(80) 



In the complex plane s the integration contour can be 
turned clockwise for 90° to be made parallel to the real 

axis. At this new integration contour s = 



s', where s' = Re(s) is a new real integration variable 
changing from to 00: 



F(t_) oc 



1 + Vw 



itt exp 



ds' e 



3/2 



t- 

TQ 



-1 + Erf 



TO 



(81) 



where, as usual, all constant and phase factors are 
dropped. Correctness of the last transformation in Eq. 
(f8"Tj) (transition to the error function) can be checked, 
e.g., by direct calculation of the integral in "Mathemat- 
ica". The normalized squared absolute value of F(tJ) 
(jSTj) is plotted in Fig. QT] as a function of the dimension- 
less variable t_/ro- As the dependence of the function 



1 / 


\ |F(t_)| 2 /|F(0)| 2 


/ 0.5 





0.5 





t-/r 



0.5 



FIG. 11: Squared absolute value of the expression on the right- 
hand side of Eq. 18H . 

F(t-) on all parameters of the crystal and pump is con- 
centrated only in the scaling factor To (|79|) . the curve of 
Fig. [11] is universal and valid for all values of these pa- 
rameters under the only condition that pump pulses are 
long, t 3> LA/c. In particular, this curve and the ana- 
lytical results of Eqs. I[77[) and (|8"Tj) arc perfectly valid for 
the above considered case of 2-ps long pulses, for which 
the temporal wave packets are described in Fig. [S] 

Found from Eqs. (|77|) and (f8lj) siglc-particle and co- 
incidence widths of the temporal distribution of photons 
are given by 



At 



(s) 



t and A< 



(c) 



0.555t 



0.555 8BLX 



(82) 



where 0.555 is the FWHM of the curve in Fig. [TTJ In 
contrast to the case of short pulses, in the regime of long 
pulses of the pump wave the coincidence width At\ does 
not depend on t 2 . For these reason, Eqs. ([82]) can be 
used for determining the entanglement parameter R t in 
the temporal representation: 



Rt 



At[ s) 



(c) 



0.555 V 8BLX0 



0.75R 



long 



0.94/f long , 
(83) 
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where i?i on g is the i?-paramctcr ([26]) calculated for long 
pulses in the frequency representation and K\ ong ([4T)]) is 
the above-calculated Schmidt number in the long-pulse 
asymptotic regime. As seen from Eq. (|53"]) , the tempo- 
ral i?-parameter, R t , has all the same functional depen- 
dences as -Riong and -Kiong and slightly differs from them 
in numerical coefficients. Moreover, R t appears to be 
strikingly close to the Schmidt number K\ ong , and this 
confirms once again a good quality of the i?-parameter 
as the entanglement quantifier. 

The coincidence width Ai^ (|82|). or equal to it "scal- 
ing factor" t ([75)) , can be interpreted as the correlation 
time of photons in the biphoton pair. It's important to 
emphasize that this time does not depend on the pump- 
pulse duration r. Hence, the correlation time remains 
finite even in the case of an infinitely long pulses or a 
purely monochromatic pump. The only reason for a fi- 
nite time of correlation in the case of long pump pulses 
is the dispersion. If the dispersion constant B formally 
tends to zero, the correlation time To tends to zero too. 

Note also that, owing to the assumption that pump 
pulses are long, in this limit the temporal walk-off ap- 
pears to be completely switched off. The scaling dura- 
tion To (|75|) . as well as the width and all structure of the 
curve in Fig. [11] arc determined only by the dispersion 
constant B, rather than by any combination of B and A 
as it was shown to occur in the case of short pulses. Fi- 
nally, the structure of the curve in Fig. [5] reminds quite 
strongly the structure of the coincidence curves for short 
pulses in the region of the peak [Fig. 0(6)] and in the 
zone III of Fig. [TTJJ This resemblance is related to the 
increasing role of dispersion and decreasing role of the 
temporal walk-off in the cases when a short-pulse pump 
exits the crystal. 



VII. CONCLUSION 

Let us summarize here our main results. 

1. For the degenerate collinear type-I SPDC process 
the parameter _R, characterizing the degree of spectral en- 
tanglement, is found and shown to be large at any values 
of pulse duration r of the pump. Even in the most unfa- 
vorable conditions r ~ 1 ps, where R(r) has a minimum, 
its value does not fall below 70, i.e., R{t) > R ni [ n » 1. 

2. The control parameter r/ (|14|) separating the regions 
of short and long pulses is found to be determined by the 
ratio of the double pump-pulse duration to the difference 
of times during which a crystal is crossed by the pump 
and idler/signal photons. 

3. By using fundamental features of the biphoton spec- 
tral wave function we managed to find the Schmidt num- 
ber K in the cases of short and long pulses (compared to 
1 ps), though the wave function was taken in the form 
significantly different from the double-Gaussian one. By 
comparing R and K, we found that they are very close 
to each other. All functional dependences of R and K 
are identical, and the difference in numerical coefficients 



does not exceed 20%. For the class of bipartite wave 
functions we consider, we found a simple universal for- 
mula [Eq. (gSJ)] establishing relation between R and K 
for any values of the control parameter rj. This means 
that as soon as the parameter R is measured and r) is 
evaluated, the Schmidt number can be easily found too. 

4. In the temporal picture, a structure of the two-time 
biphoton wave packet [^(ti, ^2) | 2 is investigated both 
numerically and analytically (the latter - with the help 
of a series of very simple formulas we have derived), in 
both cases of short and long pump pulses. It's found 
that in these two cases (short and long pump pulses) the 
shapes and localization regions of the temporal bipho- 
ton wave packet are significantly different. In the case of 
long pump pulses the localization region of [^(ti, <2)| 2 
is more or less usual: it looks like a long and thin cigar 
elongated along the diagonal t\ = t<i in the plane (t\, t%) 
(Fig. EJ6)). In contrast to this, in the case of short pump 
pulses, the localization region of |\£(ti, *2 ) | 2 is found to 
be wedge-shaped (Figs. [51 and HO]). Again, this "wedge" 
is elongated along the diagonal t\ = But, very un- 
usual, its width in the perpendicular direction (t\ = —£2) 
changes in very wide limits. In our example used above 
for estimates, from about 1.5 ps in a wide part of the 
"wedge" down to a few femtoseconds in the narrow part. 

5. The above-described structure of the temporal wave 
packet in the case of short pump pulses gives rise to the 
following specific features of the single-particle signals, 
(a) Duration of the single-particle signal is much longer 
than that of the pump (~3 ps in our example compared to 
t = 50 fs). (6) the front wing of the single-particle signal 
is anomalously long, ~ 1 ps, whereas the rear wing is as 
short as the pump pulse itself, (c) In its middle part, 
the envelope of the single-particle signal has a form of a 
plateau, i.e. is constant. The long duration of the single- 
particle signal is related to the temporal walk-off effect 
only, whereas its plateau shape and a long front wing are 
related to a combined action of the temporal walk-off and 
dispersion. 

6. In the same case of short pump pulses, coincidence 
signals have duration which depends strongly on the ob- 
servation time. Soon enough after the first signal/idler 
photons arrival to the end of the crystal, duration of coin- 
cidence signals (equal to the characteristic time of photon 
correlation in a pair) is very long, (~ 1 ps), whereas close 
to the end of observability the duration of coincidence sig- 
nals correlation time become very short (a few fs). This 
uncertainty in the definition of the coincidence duration 
makes the temporal picture in the case of short pump 
pulses very inconvenient and inappropriate for defining 
the i?-parameter and evaluating in such a way the de- 
gree of entanglement of a biphoton state. 

7. In the case of long pulses the wave packet 
*2 ) 1 2 is concentrated along the diagonal t\ = t 2 

and is very narrow and long. Duration of coincidence 
signals is small, determined purely by dispersion, and 
independent of the observation time, i.e. constant. Ow- 
ing to the last reason, the parameter R can be clearly 
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and unambiguously defined, and it appears to be very 
close to the Schmidt number. Actually, a high degree 
of closeness between the parameters R and K repeats 
itself symmetrically in the cases of (a) short pulses and 
frequency representation and (6) long pulses and tempo- 
ral representation. In both cases, with a good accuracy, 
R = K > 1. 
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